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The knowledge of the space environment is a necessity for the design 
and operation of the space vehicle, Skylab, future space schuttle and space- 
lab. Complex physical characteristics of the property of the space environ- 
ment is essential to acquire a systematic investigation to the future 
advancement of our knowledge. 

In this report, we shall present the results of our studies during 
this period of time (March 1970 - February 1974). In this study, we have 
covered a wide range of problems in the space environment, i.e., the 
problems of dynamical behavior of thermosphere, hydromagnetic wave propa- 
gation in the upper F2-region, and interplanetary space environment. Most 
of these studies have resulted in journal publications, which have previously 
been submitted to the National Aeronautics and Space Administration/Marshall 
Space Flight Center. They are as follows: 

(1) "Propagation of Hydromagnetic Waves in the Upper F2-Region," 
Planet. Space Sci. , Vol. 21, 1973. 

(2) "Kinetic Theory Analysis of Solar Wind Interaction With 
Planetary Objects," in Photon and Particle Interaction 
With Surfaces in Space (Ed. by R. J. L. Grard) , D. Reidel 
Publishing Company, Dordrecht-Holland , 1973. 

(3) "The Dynamical Responses of the Thermosphere Due to a Geo- 
magnetic Storm," AIAA Paper No. 74-217. 

(4) "An Analysis of the Upper Atmosphere Wind* Observed by 
LOGACS," Planet. Space Sci., 1974. 
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AN ANALYSIS OF THE UPPER ATMOSPHERIC WIND OBSERVED BY LOGACS 

Abstract - Wind velocities at 140-200 km altitude were observed by a 
Low-G Accelerometer Calibration System (LOGACS) flown on an Agena satellite 
during a geomagnetic storm. An interesting wind reversal observed by the 
satellite at auroral latitudes is satisfactorily explained by the neutral 
air motion caused by the E x B drift deduced from the ground-based geo- 
magnetic data recorded at stations near the meridian of the satellite orbit. 

i. INTRODUCTION ' 

High-resolution wind observations in the lower. thermosphere were made 
by an Agena satellite in an experiment designed to measure very small 
accelerations in the orbital motion. Data from this experiment designated 
LOGACS (Low-G Accelerometer Calibration System) were analyzed to determine 
the magnitude and direction of the wind component perpendicular to the 
orbital plane of the satellite (DeVries, 1972). Observations of the 
lateral wind component at 140-200 km altitude on 26 May 1967 for different 
orbits of the satellite are shown in Fig. 1. As seen for the 54th and [ 
56th orbits, a reversal of the lateral Wind direction occurred at about 
65° N latitude. In order to find a possible cause of this wind reversal, 
geomagnetic data obtained at a chain of stations near the satellite orbit 
were examined. As shown in Fig. 1, Dixon Island (DI; 73°33' N, 80°34‘ E), 
Sverdlovsk (SV; 56°44' N, 61°04- E), Tehran (TP; 35°44* N, 51°23- E) and 
Quetta (QU; 30°11‘ N, 66°57' E) are located near the meridional plane of _ 
the 56th orbit. The electromagnetic drift speeds were estimated from 
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the geomagnetic variations at these stations, and were compared with the 
LOGACS wind measurement for the 56th orbit. Since the geomagnetic D 
component during disturbed periods is strongly controlled by field- 
aligned electric currents, the H component has mainly been studied to 
estimate the electric current intensity in the ionospheres assuming that 
a significant part of this component arises from ionospheric sources. 
Whether this assumption is justifiable as a first approximation is one 
of the purposes of the present research note. 

.V.- 2. METHOD 

’ As Given by Chapman and Bartels (1940), the sheet electric current 
intensity i in the ionosphere can be estimated approximately by the 
geomagnetic variation field aB as follows: 

' . ' . I x * AB y (f/2ir) 

ly = -AB X (f/2ir) 

where f is a parameter used to eliminate the induction current effect 
(commonly assumed to be 0.6), and the subscripts x and y represent the 
southward and the eastward component, respectively. When the geomagnetic 
H component only is considered, we have 

I y — AH (f/2n) = 0.0955 AH . 


( 1 ) 
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As discussed by many scientists ( e.g . , Maeda and Kato, 1966; Matsushita, 
1967), the southward electric field E can be calculated by the follow- 
ing relation when I is ignored: 

i - • ■ - . • ■ ■ . 

E x = -I y E xy^ E xx + E yy * s xy ^ 

; v ; - -0.0955 ah z xy /(r xx i yy + r xy 2 ) (2) 

where r is the height-integrated electric conductivity. The eastward 

Ion speed V due to E x B drifts above 140 km altitude can be obtained by 

. ■ • * * . 

v„ = (E/B) cosec* :• (3) 

J' * 


where B is the earth's main magnetic field and $ is the magnetic dip 
angle. 

The neutral wind velocity V n can be approximately determined from 
the ion drift velocity V. as follows; 

1 n 








where |P n is the stress tensor, p n is the neutral density, and T pi is the 
neutral-ion collisional time. If we neglect viscosity and assume that the 
neutral pressure gradient can be expressed in terms of neutral velocity as, 




then (4) becomes 



dV 

n _ 




with-a = a and g being constants indicating the order of magnitude 

of the pressure gradient. The solution of (6) is: 



v i (f) 




This expression indicates that the neutral wind velocity approaches the 
ion drift velocity with a delay time of approximately T m -/0“ a )» which 
is on the order of 5x10^ sec for |a| « 1 ( e.g . , Banks, 1966; Fedder and 

^ ^y -y 

Banks, 1972). In other words, through ion-neutral collisions, V n - V.. = 

may occur within a few hours. When V lasts for a few hours, then, a neutral 

«/ 

wind flow which has approximately the same velocity as may be generated. 
Accordingly, V estimated from AH using Equations (1), (2) and (3) can 
be compared with the LOGACS neutral -wind measurements, taking into con- 
sideration a few hours time lag. 




As presented in Fig. 2, DI shows a positive variation in the H Fig. 2 

component, while the three other stations (SV, TP and QU) in lower 
latitudes than DI show a negative variation, at about 6h UT on 26 May 
.1967 corresponding to the time of the 56th satellite orbit. To estimate 
AH, the zero level is assumed to be the average value of 6-7h UT on 
25 May as shown by the horizontal straight line in Fig. 2, Since the aH 
due solely to ionospheric electric currents is needed to discuss ionospheric 
' winds, ring-current effects must be eliminated. From the equatorial d 
values obtained by Sugiura and Cain (1970), ADst«cose was computed 
.. (broken curves in Fig. 2), where e is the latitude angle of each station 
and AD st is the deviation from the d t value averaged' for 6-7h UT on 
25 May. ♦ " / 

Deviations of H values from the zero level (solid circles) and those 
from ad '* cose (crosses) for the four stations at 6h UT on 26 May are 
presented in the top diagram of Fig. 3. Solid circles indicate the Fig. 3 

extreme case with no Dst contribution. It can be noticed that the sign 
of AH changes at 60°-70° N latitude, roughly corresponding to the change- 
over location of the LOGACS wind direction shown in Fig. 1. 

For simplicity, the' scaling of AH is made for 6h UT on 26 May correspond- 
ing to the 56th orbit. As discussed at the end of the preceding section, a 
few hours are required to produce a neutral wind velocity of the same order 
as the ion drift velocity. Since the aH values remain generally the same 
during a few hours before 6h UT on 26 May (see Fig. 2), aH at 6h UT is , 
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acceptable as a representative value of the geomagnetic variation for the 
present discussion. 

To obtain E x from Equation (2), values of the height-integrated 
electric conductivity E are necessary. Since these values during geo- 
magnetic storms are difficult to estimate accurately, the S values during 
quiet periods (Tarpley, 1969) are substituted, bearing in mind that the 
resulting |E I mqst be close to the maximum value. The adopted con- 
ductivities at four geomagnetic stations are listed in Table 1. Cal- |_Ta_ 
culated E x distributions with (crosses) and without (circles) Dst effect 

are shown in the middle of Fig. 3. 

With E x from Equation (3), the eastward drift speed V y can easily 

be obtained. These V y are compared with the observed neutral wind speed 
by LOGACS in the bottom diagram of Fig. 3. Agreement between the two is 
satisfactory as a first order approximation. However, the absolute values 
of the LOGACS wind speeds at latitudes higher than 50° N are slightly 
larger than the calculated |V y | values which are close to the maximum 
since they are obtained from the near maximum ]E X | values. One possible 
explanation of this discrepancy is that |aH| recorded at the ground is 
smaller than that which can be attributed to an electric current of purely 
ionospheric origin. Field-aligned electric currents and some magnetosphenc 
variations may cause a reduction of 1 aH| at the ground. Although various 
data and careful examination are needed to discuss the details, an 
interesting conclusion is that the reversal of the LOGACS wind direction 
near 65° N is understandable in terms of the reversed ionospheric drift 
directions between Dixon Island and Sverdlovsk inferred simply from the. 
geomagnetic H records at these stations. 
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... TABLE 1. 

HEIGHT- 

INTEGRATED ELECTRIC CONDUCTIVITIES 


’ 

(in 10" 9 emu) 


* 

’ 2 xx 

£ xy 

s yy 

DI * 

' 4.27 

. 5.38 

• 

4.19 

SV 

5.79 

7.42 

5.24 

TP * 

11 .5 

1*4.0 

78.3 

*r-.. 

QU 

' 15,5- 

17.8 

89.0 
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FIG. 2 


FIG. 


Figure Captions 

. Distributions of wind velocities normal to the satellite path at 
about 150 km altitude measured by LOGACS for its 53rd, 54th and 
56th orbits on 26 May 1967 are shown in geographic coordinates 
(with three Greenwich longitudes and the equatorial semi-circle). 
The 56th orbit was along the 60°-70° E meridian at approximately 
6h UT. locations of four geomagnetic stations (Dixon Island, 
Sverdlovsk, Tehran and Quetta) near the 56th orbit are shown by 
solid circles. ; 

. Geomagnetic H variations (positive upward) recorded at four 
stations on 25-26 May 1967. The zero level is assumed to be the 
average value of 6-7h UT on 25 May, and shown by the horizontal 
straight line. The broken curves indicate ADst*cose measured 
from the zero level, whereADst is the deviation from the Dst 
value (+30y) averaged for 6-7h UT on 25 May and 9 is the latitude 
angle of each station. 

. ’• t 

n 

, (top) Deviations of H values from the zero level (solid circles) 
and those from ADst«cose (crosses) at 6h UT on 26 May 1967. 
(middle) Southward electric field, E^, calculated from H 
deviations. 

(bottom) Observed neutral wind speed by LOGACS (small solid 

circles with a fitted smooth curve obtained from the wind speeds 
• 

for the 56th orbit in Fig. 1) are compared with calculated east- 
ward drift speed V (large solid circles and crosses). 










